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Abstract 

The nine-component positive vector optical tomographic probability portrait of quantum state of 
spin-1 particles containing full spatial and spin information about the state without redundancy is 
constructed. Also the suggested approach is expanded to symplectic tomography representation and 
to representations with quasidistributions. The evolution equations for constructed vector optical and 
symplectic tomograms and vector quasidistributions for arbitrary Hamiltonian are found. The evolu¬ 
tion equations are also obtained in special case of the quantum system of charged spin-1 particle in 
arbitrary electro-magnetic field, which are analogs of non-relativistic Proca equation in appropriate 
representations. The generalization of proposed approach to the cases of arbitrary spin is discussed. 

The possibility of formulation of quantum mechanics of the systems with spins in terms of joint prob¬ 
ability distributions without the use of wave functions or density matrixes is explicitly demonstrated. 

Keywords: Quantum tomography, spin tomography, evolution equation, Proca equation, non-negative 
vector portrait of state. 
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1 Introduction 

The proposition of optical tomographic description of states of spinless quantum systems was formulated 
in mm- Generalising the optical tomography technique the symplectic tomography was suggested, and 
evolution equation for symplectic tomograms of spinless quantum systems was found in mm providing 
a bridge between classical and quantum worlds (for review see [5]). Evolution equations for optical 
tomograms of spinless quantum systems were obtained in Refs. mm- 

In Ref. [8] the spin density matrix for particles of arbitrary intrinsic angular momentum is explicitly 
expressed in terms of directly measurable expectation values of components of multipole moments, or 
the relative weights of partial beams split-up by a Stern-Gerlach apparatus. 

Extending the tomographic approach in Ref. [9] the spin tomography was formulated based on the 
description of spin states with the help of positive distribution functions depending on continuous variables 
like Euler’s angles, a spin state reconstruction procedure similar to the symplectic tomography was 
considered, and quantum evolution equation of spin dynamics was found for continuous spin tomogram. 

In Ref. [To] spin dynamics was expressed for expectation values of spin projections along a discrete 
set of fixed directions. The spin tomography was also studied in [miiaiiiKiiKisKinKiTKisKini and in 
other papers. 

The first attempt of tomographic formulation of the Pauli equation simultaneously describing both 
spatial and spin dynamics, apparently, was done in [20]. The evolution equation obtained is extremely 
complicated, because it uses redundant tomogram depending on continuous Euler’s angles and on sym¬ 
plectic vareables. 

In our resent paper m we introduced the positive vector optical tomogram fully describing both 
spatial and spin characteristics of the quantum state of spin-1/2 particle without any redundancy. We 
obtained the evolution equation for this vector optical tomogram and considered examples of evolution 
of quantum systems in proposed representation. Also we discussed the expansion of our approach to 
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representations of Wigner and Husimi quasidistributions (and pointed out the possibility of dissemination 
of the discussed scheme to the Glauber-Sudarshan representation). 

The aim of our work is the construction of spin-1 particle quantum state vector tomography without 
redundancy of information representing the joint vector distribution for space coordinates and spin pro¬ 
jections; and derivation of the evolution equation for such distribution, which would be an analogue of 
non-relativistic Proca equation. The latter will allow to explicitly demonstrate the possibility of formu¬ 
lation of quantum mechanics of the systems with spins in terms of joint probability distributions without 
application of wave functions or density matrixes. 

The paper is organized as follows. In Sec. [2] we give basic formulae of tomographic and quasiprobability 
representations of quantum mechanics for spinless particles. In Sec. [3] we introduce a positive nine- 
component vector probability and quasiprobability description of spin-1 particles and give the evolution 
equations for such vector-portraits of quantum state with arbitrary Hamiltonian. In Sec. [4] charged spin- 
1 in arbitrary electro-magnetic field is considered in proposed representations, and evolution equations, 
which are analogs of nonrelativistic limit of Proca equation, are obtained. The conclusion is presented in 
Sec. [5l 


2 Probability representation and evolution 
of spinless quantum systems 

Let us review the constructions of tomographyc or quasidistribution representations in general case for 
spinless systems. 

If the state of the quantum system is described by the density matrix p normalized by the condition 
Trp = 1, then in accordance with general scheme the tomographic distribution function or quasidistribu¬ 
tion w{x,r]) is related with the density matrix as follows (see [22]): 


w{x,p,t) = Ti{p{t)U{x,p)}, pit) = j wix,r},t)Dix,p)dx dp, 


( 1 ) 


where x is a set of distribution (quasidistribution) variables, p is a set of parameters of corresponding 
tomography, and Uix,p), Dix,p) are dequantizer and quantizer operators for appropriate tomographic 
scheme or corresponding quasidistribution representation. For Wigner [23|, Husimi [23], and Glauber- 
Sudarshan |25[ l26] quasidistributions the set of parameters p is absent. 

Notion of quantizer and dequantizer is related to star product quantization schemes (see recent review 

m)- 

Quantizer and dequantizer are constrained by the duality relation 

Tr{?7(x, p)D{x', p')} = 6{x — x')6{p — p'). 

The von-Neumann equation without interaction with the environment 


d 


(2) 


with the help of maps of type ([T]) transforms to evolution equations for tomograms [7], or to Moyal 
equation [28] for the Wigner function [23], or to evolution equation for other quasidistribution 


dtw{x,p,t) = — / Im 

n I 


Tr 


^H{t)Dix',p')Uix,p)'^ w{x',p',t)dx'dp'. 


( 3 ) 


If we have spinless quantum system in the V—dimensional space, then dequantizer and quantizer for 
optical tomography equal 


N 


u^{x,9) = \x,e){x,e\ = \[5[x„-q ^cos 


cr=l 


sin 6„ 




( 4 ) 
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N 


bw{X,9)= I - q^cos9a - \ d^r], 


a=l 


2TTm„uJn 


rriaUJa 


( 5 ) 


where nia and cOa are constants that have the dimensions of mass and frequency and are chosen for reasons 
of convenience for the Hamiltonian of a quantum system under study, \X,9) [22] is an eigenfunction of 

the operator X{9) with components cos 9^ + {p^ sin 9^^) /(mo-cuo-) corresponding to the eigenvalue 

X. 

For symplectic tomography one can be written 


N 


{,X\X^ ) (X, /i, I (5(Xfj quPfj Pa^u '); 


( 6 ) 


(J=l 


N 


<j=l 


bM{X, fl,u)=Y\ exp <1 i 


nriaUJa 


h 


(Xq- qcrl^cr f ? 


( 7 ) 


where |X,/I, F) is an eigenfunction of the operator X{jl,v) with components X^ = PaQa + r'o-Po- corre¬ 
sponding to the eigenvalue X. 

For Wigner representation we have 


Uw{(i,p) = 


1 


|q — u/2) exp(—ipu/h)(q -|- u/2| d'^u, 


(IttH) 


N 


( 8 ) 


p) ^ y exp(2zpu/h)|q u)(q - u| ; 
for Husimi representation (see [221 [5U] 1 

1 ( 


UQ{(i,p) = {2TTfi) ^|a)(a|, a = 


mu) i 

q+ — p I, 


V2\\ h y/h'. 


muj 


where |q) is an eigenvalue of the position operator, |a) is a standart boson coherent state, 

-DQ(q,p) = (^) ^ y'd^xd^y||x)(y|exp (^^(x-y)2) 


(9) 


( 10 ) 


X exp 


N 


n 


mui / X -|- y 
h 




moj, n 2 * / N 

— (x-y)^ + -p(x-y) 


( 7=1 L?if7=o 


n, 


d2^ 


muj muj 


(11) 


Likewise, the Glauber-Sudarshan P-function 
corresponding dequantizer and quantizer 


(see, also m) can be introduced with the help of 


el“l^ 


Up{a) = ^ / l/3)(-/3|el/'l P , bp{a) = |a)(a| , 


TT' 


( 12 ) 


and so on for the other tomographic schemes. 
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3 Probability description of spin-1 particles 


In general case the evolution of charged spin-1 particle in the external electro-magnetic field is deter¬ 
mined by Proca equation [321133| . This is a relativistic wave equation of four-component wave function 
{ipo, ifi, ip 2 -, ips)- But in the case of weak relativism it can be reduced to the Schrodinger type equation 
for three-component spinor wave function with the Hermitian Hamiltonian ('i/ii,V’a) [341135] , 

For quantum system of charged spin-1 particles without electrical quadrupole moment, with charge e 
and mass m in the electro-magnetic field with vector and scalar potentials A(q, t), p{(i, t) this Hamiltonian 
has the form: 

H = — {p--A) + eif --sH = Ho --sH, 13 

where Hq is an independent on spin part of Hamiltonian, H = rotA is a magnetic field, and x is a 
magnetic moment of the particle. 

The wave function satisfy by the normalization condition 

j (IV'i(q)P + IV'2(q)P + |^/’3(q)P)d^g = i- 

So, mixed states are described by the density matrix pij with dimension 3x3, which is actually defined 
by nine real scalar components. 

Analogously with the case of spin-1/2 particles [21j . to construct the vector portrait of such density 
matrix we must solve the state reconstruction problem, i.e., we have to find the inverse map, which 
transforms the set of expectation values of observables constituting a quorum to the density matrix. 

For this purpose we should choose nine spin-1 states luj, n*) with definite spin projections cjj along the 

directions rij, which define nine-component dequantizer vector ZY of 3 x 3 spin matrices with components 
Ui = |(Tj, nj)(crj, nj| and quantizer 3x3 matrix H of nine-component vectors so, that 

2S-I-1 (2s-|-l)2 

Tr/cZ {^j{ki)T^(ki)j'} = ^ ^j{ki)'^{ki)j' = ^jj'i ^ ^j{ki)^{k'i')j = ^kk'^w- (14) 

k,l=l i=l 

Here the first index j = 1, 2,..., 9 is the number of the component of the nine-component vector, and {kl) 
are the indexes of 3 x 3 matrices. It is obvious that the set of matrixes {Ui} must be linearly independent. 

With the help of U the nine-component vector tomogram or quasidistribution is defined as 

t) = Tr |/5(t)[/(a:,r/) (8 )Zy| , (15) 

where the trace is calculated also over spin indexes, and U{x,rj) is defined by formula dH), ([6]), ([8]), 
m , or m- Here w{x,rj,t) is the aggregate designation of the vector optical w{X,9,t) or symplectic 
M(A,/2, F, t) tomogram, vector Wigner function lF(q, p, t), vector Husimi function (5(q, p,t), or vector 
Glauber-Sudarshan function P(a,t)- For optical and symplectic vector tomograms and for the Husimi 
vector quasidistribution each of the function Wj{x, tj, t) is the probability distribution of the operator x(r/) 
at time t under the condition that the particle has the corresponding value of spin projection along the 
appropriate direction. Consequently, the components of the vector w^x, rj, t) must be integrable over da; 
and must satisfy the inequalities 

0 < Wj{x,r],t) < 1, 0< J Wj{x,r],t)dx < 1, j = 1,..., 9. 

The components of the constructed vector Wigner function and vector Glauber-Sudarshan function cor¬ 
responding definite spin projection along the appropriate direction are not obligatory non-negative, but 
definition (|15ll guarantee that they are definitely real. 
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Such a definition (I15jl for our vector Wigner and Husimi functions differs from those usually given 
in literature by many authors, when the Wigner function and Husimi function Qjfc(q, p,t) 

become (2s+l) x (2s+l) matrices dependent on position and momentum, but their non-diagonal elements 
over the spin indexes are not surely real. So, the main advantage of such quasidistributions with respect 
to density matrix disappears. 

The inverse map of (I15p is written by means of spin quantizer D as follows 

Pjk{t) = / D{x,r]) ^'D(^jk)w{x,r],t)dxdr]. (16) 


Generalizing equation ([3]) to the case of spin particles we can write the evolution equation for the 
components of the tomogram or vector quasidistribution 


(2s-hl)2 


^ V ^ 

dtWj{x,r],t) = ^ 2 ^ 


Im 


k=l 


Tr |f7(x, r]) (g) UjH b{x', f/') ® 


xwk{x', r] , t)dx dr] , j = 1, (2s -|- 1)^. 


(17) 


Let us point out that the scheme proposed admits a generalization to the case of arbitrary spin s. If 
we have the evolution equation of the quantum system with spin s for the nonnegative, hermitian, and 
normalized density matrix, then we can introduce (2s-|-l)^-component vector of (2s-|-l) x (2s-|-l) matrices 

dequantizer U and dual (2s -t- 1) x (2s -|- 1) matrix of (2s -|- l)^-component vectors V, which are related 
by conditions ()14p . After that we can define (2s-|- l)^-component vector tomogram (or quasidistribution) 
in accordance with analog of (jl5h . and with the help of the formula analogous to (I16p we can write the 
evolution equation of type (|17ll for the (2s -|- l)^-component vector tomogram or quasidistribution. 


4 Example of vector tomography representation for spin-1 particle 

To determine the dequantizer we should choose nine positive projections of the quantum state, which 
completely define the density matrix. Let us choose such spin projectors as follows: 

U = Qsz = l)(s^ = 1|, |s^ = 0 )(s 2 = 0|, Is^ =-l)(s 2 =-1|, 

I'Sa; — 1)(S3; — 1|) — 0)('Sa; — 0|) \^xy — l)('Sa;j/ — 1|) 

\^xy — 0)('Sa;j/ — 0|, \Syz — 0) (Syj; — 0|, \Sxz — 0) ('Sxz — b| ^ , (18) 

where |sj = ±1, 0) is an eigenfunction of the projection of spin operator to the direction j corresponding 
to the eigenvalue ±1 or 0, and \sxy), \syz)-, are eigenfunctions of projections of spin operator to 

directions exy = 0), eyz = (0, l/\/2, l/\/2), &xz = (l/\/2, 0, l/\/2) respectively. 

Choose the spin representation, in which components of spin operator are defined as follows: 



' 0 

1 

0 ' 


' 0 -1 

0 


' 1 

0 

0 

1 




i 








1 

0 

1 

’ ""^71 

1 0 

-1 

II 

0 

0 

0 

_ 0 

1 

0 _ 

. 0 1 

0 


_ 0 

0 

-1 
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After calculations in matrix notations for dequantizer lA we have 


U = 


1 0 0 
0 0 0 
0 0 0 


0 0 0 
0 10 
0 0 0 



' 0 

0 

0 ' 

1 

’ 4 

1 

72 

1 


0 

0 

0 

72 

2 

72 


_ 0 

0 

1 _ 


1 

72 

1 


1 

1 0 -1 

1 

1 

1 — i —i 

1 

1 0 z 

2 

0 0 0 

-1 0 1 

’ 4 

z + 1 

i 

2 1 -z 

l + i 1 

’ 2 

0 0 0 

-i 0 1 


1 

4 


1 1 
-iy/2 2 —i-\/2 

1 iV2 1 


1 

’ 4 


1 -V2 -1 

-V2 2 V2 

-1 1 


(19) 


From duality relation (jl4p after some calculations we obtain spin quantizer D, which is a 3 x 3 matrix of 
nine-component vectors 


( 20 ) 


where ( jk ) are the indexes of 3 x 3 matrix and I = 1,2, ...,9 is the index of the component of nine- 
component vector 


- r - 1 

^(11) 

7(12) 

7(13) 

^ = {^(ifc)«} = 

7(21) 

(^(22) 

(^(23) 


- ^(31) 

(^(32) 

^(33) 


V 


( 11 ) 


V 


( 12 ) 


V 


( 13 ) 


V , 


( 22 ) 


V 


( 23 ) 


V 


( 33 ) 


1, 0, 0, 0, 0, 0, 0, 0, o), 


1 I-V 2 I-V 2 1 I-V 2 1 + i 

~VT' 

L\ 


l+i . i i 

^’"*’" 2 ’7!’ 72 ^’ 

/I — i I — i \ 

0, 0, -1, 0, z, 0, oj, 

0, 1, 0, 0, 0, 0, 0, 0, 0), 

1 i 1 i 1 zl-|-z , i i 1\ 

"^^2’ "71 + 2’ "771 + 2’ "2’ "71’ 71J’ 


(^ 0 , 0 , 1 , 0 , 0 , 0 , 0 , 0,0 


^( 21 ) - ^ 7 )’ ^( 31 ) “ ^ 7 )’ ^( 32 ) “ ^* 23 )- 

Obviously that for such of definition (|18p of dequantizer U, three components of the vector w{x,rj,t) are 
normalized by the condition 


wi{x,ri,t)dx + J W 2 ix,r],t)dx + J W 3 {x,r],t)d. 


X = 1. 


( 21 ) 


For Hamiltonian (1131) the evolution equation (11711 of the optical vector tomogram is written as follows 
(see m): 

dtw { x , e , t ) = MUx , e , t ) w { x , e , t ) + sux , e , t ) w { x , e , t ), (22) 
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where 


= ^ImHo ([q]^(X,0),[p]^(X,0),i 

is an operator depending on position [q]^ and momentum [p]^ operators in the optical tomographic 
representation 


[q^]yj{X,e) = 


\pa\u,{X,9) = mOJa -cos 00 


d 


n -1 


dX„ 


+ Xfj cos Oa + i 


. h sin Ofy d 


2mfjL0o- dXa ’ 


d 


dXa 


-1 


d ^ ■ n \ /I ^ 

^ + --cose„—. 


(23) 


(24) 


and Stp(X, 0, t) is a 9 x 9 matrix operator, responsible for the interaction of spin with the magnetic field 

2s+l 


2^ 


sH [q]^(X,0 


T^{m'l)k 

J {mm') 


(25) 


With omitted arguments and introduced designations 

[i,]^ = ( [qUX, 0),t), H, = [H,]^ = H,- ([q]^(X, 0 ), t) , 

VqAl = VqA(q^ [q]^(X,0),t) 

the explicit form of A4^ in general case of time-dependent and non-homogeneous electromagnetic field is 
written as 

3 


M^{x,e,t) = ^ 


a;ri 


n=l 


COS 


3 1 f 3 

^ ^ — - sin 20n < 1 -I- 




’"d0n 2 ”1 

, r 1 19 2e ^ r I -.1 e 


2 e 

h 


Im [(^]t 


■Im[A]^--Im Ap H-Re[VqA] 

mc^h men L Jiij me 


For symplectic vector tomography we can find the evolution equation 

dtM{X,fI,n,t) = MMiX,fI,u,t)M{X,fI,n,t) + Sm(A,/ 2, F, t) M(A,/T, F, t), 

where operator A4 m(A,/T, F, t) corresponds to spinless part Hq of the Hamiltonian (fTH)) 

2 ~ ^ \ d 2e 

MM{X,fl,u,t) = -lmHo{[p]MiX,fI,u),[q]MiXfI,u),t) = fl— +—lm[ip]M 


(26) 

(27) 


+ 


^ T rAi 2 2 e 


Ap + —Re[VqA]j^^, 

J M me 


(28) 


where 


[Aj]M = Aj [[c^M{X,fl,v),t) , 

[VqAjjVf = VqA ^q - 

and [q]M, [p]m are position and momentum operators (1^^ in the symplectic representation (see m) 


\M = ^[ [qjM(A,/2,z7),t) , 
]MiX,il,n),t) , 


\Pa 


\M = 


d 


dX, 


-1 


d ./Xo-h d 

durr * 2 dX„ 


[i 


a\M — 


d 


dX, 


-1 


5/i, 


d u„h d 
+ 1 - 


2 dX, 


(29) 
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The 9x9 matrix operator is defined by the similar formula (I25p . where the operators of 

components of the magnetic field Hj must be replaced with corresponding operators in the symplectic 
tomography representation [HjJm = Hj . 

Making similar calculation we can obtain such evolution equation for our vector Wigner function, 
which is a generalization of the Moyal equation [28] 


A 

¥t 


iy(q,p,t) = 


p (9 2e 


+ - 


m (9q h 
2 e 


, ih d \ e 


, ih d 

|q+y^,‘ 


mch 


Im^A|q+|A,t 


ih d 


+ 


me 


ReVqA q —> q + 


ih d 
2 9p’ 


,t] + Siv(q,p,t) 




(30) 


where 9x9 matrix operator Siy(q)P)i) is defined by the same formula (|25l) . where the operators of 
components of the magnetic field Hj must be replaced with corresponding operators in the Wigner 


representation H, (q + 


ih _d_ 
2 dp 




The corresponding generalization of the evolution equation of the Husimi function [36| to the case 
of vector quasidistribution has the form (for simplicity we choose the system of measurements so that 
m = oj = h = 1): 


m 




+ 


d I d d 2e ^ ( Id id 

- 2 8^5 + X+ 2 85 + 2 85'* 

^ '' 1 d i d 


^ ImA^ ( q + 


2 5q 2dp^ 


2e 


Id id 


- -'"■i^l‘>+285+285’‘ 


15 id 
2^ “ 2^ 


+ 


ReVqA I q —> q 


15 id 

2 5q ^ 2 5p ’ 


,t] + SQ(q,p,t) 


Q(q,p,t), 


(31) 


where 9x9 matrix operator SQ(q, p, t) is defined by (1251) in which components of the magnetic field Hj 
are replaced with Hi(q+^^ + §^ 


,t . 


5 Conclusion 

To resume we point out the main results of our paper. We have constructed the nine-component positive 
vector optical tomographic probability portrait of quantum state of spin-1 particles, which contains full 
spatial and spin information about state without redundancy. We have expanded suggested approach to 
symplectic tomography representation and to representations with quasidistributions. All of the compo¬ 
nents of the constructed vector Wigner function are real, and all of the components of the vector Husimi 
function are non-negative. 

We found the evolution equations for such vector optical and symplectic tomograms and vector 
quasidistributions for arbitrary Hamiltonian and obtained these equations in special case of the quantum 
system of charged spin-1 particle in arbitrary electro-magnetic field, which are analogs of non-relativistic 
Proca equation in appropriate representations. Also we discussed the generalization of our approach to 
the cases of arbitrary spin. 








The general equations obtained are relatively complicated, but in many special cases they are much 
simpler and could allow for the possibility of analytical and numerical solutions. 

The results of the paper explicitly demonstrate the possibility of formulation of quantum mechanics 
of the systems with spins in terms of joint probability distributions without application of wave functions 
or density matrixes. 
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